Distribution of Fluctuational Paths in Noise-Driven Systems* 

M.I. Dykman and V.N. Smelyanskiy 
Department of Physics and Astronomy, 
Michigan State University, East Lansing, MI 48824, USA 

Abstract 

Dynamics of a system that performs a large fluctuation to a given state is essentially determin- 
istic: the distribution of fluctuational paths peaks sharply at a certain optimal path along which 
the system is most likely to move. For the general case of a system driven by colored Gaussian 
noise, we provide a formulation of the variational problem for optimal paths. We also consider the 
prehistory problem, which makes it possible to analyze the shape of the distribution of fluctua- 
tional paths that arrive at a given state. We obtain, and solve in the limiting case, a set of linear 
equations for the characteristic width of this distribution. 



* This paper is dedicated to Rolf Landauer, with deep respect and admiration 
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I. INTRODUCTION 



Large fluctuations, although infrequent, play a fundamental role in a broad range of 
processes, from diffusion in crystals to nucleation at phase transitions, mutations in DNA 
sequences, and failures of electronic devices. In many cases the fluctuating systems of 
interest are far from thermal equilibrium. Examples include lasers, pattern forming systems 

trapped electrons which display bistability and switching in a strong periodic field 0, [J, 
and spatially periodic systems (ratchets) which display a unidirectional current when driven 
away from thermal equilibrium |4|. 

nn 

It was very clearly shown by Landauer [5|, [6( that, whereas for systems in thermal equi- 
librium the probabilities of fluctuations are known at least in principle, for nonequilibrium 
systems there is no universal relation from which these probabilities can be obtained: even 
though the system mostly stays in the vicinity of one of locally stable states, the distribution 
over these stable states can be found only from global analysis. This distribution may be 
strongly affected by nonthermal perturbations in the rarely occupied intermediate states, 
i.e. by the large fluctuations which determine the probabilities of switching between the 
stable states (Landauer's blowtorch theorem). 

The major physical problems in the theory of large fluctuations are not only calcula- 
tion of the fluctuation probabilities, but also analysis of the dynamics of large fluctuations. 
Understanding this dynamics is particularly important for controlling large fluctuations. 

An intuitive approach to the theory of large fluctuations makes use of the optimal path 
concept. The optimal path is the path along which the system is most likely to move when 
it fluctuates to a given state from the vicinity of the stable state. Although trajectories 
of a fluctuating system are random, it is clear from Fig. 1 that the probabilities for the 
system which is found at a point qj at an instant tf to have arrived to this state along 
different paths are very different: e.g., it is unlikely that the system has been staying far 
from the equilibrium position for a long time, or that it has experienced an extremely large 
acceleration, as is the case for the paths 1 and 3 in Fig. 1. 

Optimal paths and the probability distribution of the fluctuational paths are real physical 
objects: they have been observed in analog experiments ?| ls[ I^] and digital simulations 10]. 
In the theory of large fluctuations, the pattern of optimal paths plays a role similar to that 
of the phase portrait in nonlinear dynamics. 
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The fundamental role of the distribution of fiuctuational paths was recognized already 



by Onsager and Machlup 



in fact, they obtained optimal paths for a linear system 



in thermal equilibrium with the bath with a short correlation time (the approximation of 
Brownian motion). For such systems, whether they are linear or nonlinear, the optimal path 
to a given state is the time-reversed path from this state to the vicinity of the stable state 
in the neglect of fluctuations (the deterministic path) [9], [l2]. This is no longer true for 
nonequilibrium systems, because, in general, they lack time reversibility. Even 
nonequilibrium systems the pattern of optimal paths may have singular features 



or simple 

3,Q. 



U(q) 





FIG. 1: (a) A particle with a coordinate q fluctuating away from the bottom q = q^ of the potential 
well U(q) to a point qf. (b) Various paths along which a particle can reach qj at a given instant 
tf in the course of the fluctuation; the probability densities for moving along different paths are 
exponentially different. With overwhelming probability the system moves to a point (qf,tj) along 
an optimal path. 



In the present paper we provide a general formulation of the problem of optimal paths 
for nonequilibrium systems driven by Gaussian noise. In Sec. II, based on a path-integral 
expression for the transition probability density, which allows for a prefactor, we derive 
an integral variational functional for optimal fiuctuational paths. In Sec. Ill we provide a 
formulation of the prehistory problem for the distribution of fiuctuational paths to a given 
state. This formulation is reduced to a linear integro-differential equation; the explicit form 
of this equation depends on the form of the correlation function of the noise. In Sec. IV we 
discuss the solution of this equation in the case where the destination state is close to the 
stable state of the system. Sec. V contains concluding remarks. 
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II. LARGE FLUCTUATIONS INDUCED BY GAUSSIAN NOISE WITH AN AR- 
BITRARY POWER SPECTRUM 



Nonlinear systems driven by nonthermal Gaussian noise form an important and fairly 
general class of nonequilibrium systems. In spite of seeming simplicity, they display a variety 
of interesting effects, unidirectional current in periodic structures being an example [4]. 
Substantial progress in the theory of large fluctuations in such systems has been made 
within the last two decades, particularly by implementing the path integral technique (see 
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0, Q, 



19| and references therein). So far this technique has been applied to 
systems driven by the noise which is a component of a Markov process [3], in which case 
the reciprocal power spectrum of the noise l/$(u;) is a polynomial in u 2 , where 

$H= / dte^t), <f>(t) = </(t)/(0)). (1) 



It is advantageous to express the characteristics of fluctuations in the system in terms of the 
noise power spectrum, since can often be measured in experiment. The shape of <3>(u;) 
depends on the source of the noise and the coupling between the dynamical system and this 
source. 

The noise with of the form of a polynomial in u 2 , although interesting, is not the 

most general type of noise. An important example of Gaussian noise with a non-polynomial 
$ _1 (u;) is the noise with Gaussian power spectrum in the central part 

$0) = £>exp [- (lu 2 - tt 2 ) 2 /4tt 2 a 2 ] , \u 2 - Qg| < a 2 . (2) 

In particular, noise of this sort is produced by the electric field of inhomogeneously broadened 
(in particular, Doppler-broadened) radiation. 



A. Transition probability in a system driven by Gaussian noise 

In this section we provide a general formulation of the problem of large fluctuations 
induced by Gaussian noise. We consider stationary systems or systems in a time-periodic 
field. To simplify notations, we will assume that the system under consideration is described 
by one dynamical variable, q. The Langevin equation of motion is then of the form: 



q = K(q;t) + f(t), K(q; t + T) — K(q; t), 
4 



(3) 



where f(t) is the zero- mean stationary noise. We assume that the noise is characterized by 
a certain correlation time t corr over which its correlation function decays (at least exponen- 
tially, in the limit of large time). 

For weak noise intensities, over a time t which exceeds t corr and the characteristic relax- 
ation time in the absence of noise t re i, the system will approach the stable state q^(t) and 
will then perform small fluctuations about it. In periodically driven systems the state q^°\t) 
is also periodic and is given by the equation 

= K(q® ; t) , g (0) (t + T) = q® (t) . (4) 

(we assume that the period of the state q^ is the same as that of the periodic driving). 

In the course of a large fluctuation, the dynamical system is brought from the attractor 
to a distant point qf at the instant tf (cf. Fig. 1). For this to happen the system should 
have been subjected to finite forcing over certain time. Different realizations of the force 
f(t) can result in the same final state. The system trajectories q(t) for each realization of 
f(t) are deterministic rather than random, and they are independent of the characteristic 
noise intensity 

D = max $(cl>). (5) 
The probability density of realizations of f(t) is given by the functional (cf. [2^) 



V[f(t)} = exp 



~ J dtdt' f(t)F(t-t')f(t') 



(6) 



where the operator !F(t) is related to the correlation function of the noise <p(t) by the 
expression 

J dt x F(t - - f) = D5{t - If). (7) 

In some cases (in particular, for the noise f(t) being a component of a Markov process) a 
formal solution of this equation can be written as 

T(t) = D5(t)/®(-id/dt). (8) 

Here, we have taken into account that the noise correlation function is even, <p{t) = <p(—t), 
as is also the noise power spectrum, <3>(a;) = $(— uj). 

One can write the probability density p(qf,tf) for the noise-driven system to arrive at 
the point qf at the instant tf, provided it has started from the point qi at the initial instant 



U, as a path integral 



P(Qf,tj 



Vq(t)S[q(t)-q det (t;f\ qi )} 



(9) 



Here, qdet{t] f\li) is the solution of the dynamical equation of motion ([3]) for a given real- 
ization of the noise f{t), and S[q(t) — qdet{t)} is the functional delta-function: it peaks at 
the function q(t) equal to qdet{t)- The averaging (...) means integration over f(t) with the 
probability density functional |6| as a weighting factor 20j. In what follows we assume 
that ti — > — oo and the initial point is close to the attractor, qi « Q^{U)- I n this case the 
function p ([9]) gives the stationary probability distribution, which is periodic in time for the 
time-periodic force K in ([3]), p(qj,tj) = p(qf,tf + T). 



It is convenient to perform averaging over / in (Q by writing the 5-function in the 
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a path integral over an auxiliary variable k(t)/D. Using standard transformations 
one can show that the expression for p(qf,tf) can be written in the following form: 



orm of 

3,Q 



p(Qf,tf) = C J Vf(t)V[f(t)] J V^Vq(t) 

At) r . 



(10) 



x exp 



dt 



[q-K(q;t)-f(t)]--K' 



K'(q;t) 



dK{q-t) 
dq 



where C is the normalization constant, and V[f(t)] is the probability density functional for 
the random force ([6]). 



B. Variational problem for optimal paths 

If D is sufficiently small, as we assume, then, for all f(t) which result in a large fluctuation 
to a given state, the values of the probability density functional ([6]) are exponentially small; 
they also exponentially strongly differ from each other for different f{t). Thus one would 
expect that there exists one realization f(t) = f op t{t) which is exponentially more probable 
than the others. This realization provides the maximum to the functional ([6]) subject to 
the constraint that the system described by Eq. ([3]) is driven to a designated state qf. 
Respectively, there are in fact two interrelated through ([3]) optimal paths: that of the 
system, q op t{t), and that of the force, f op t{t) [15, 



Formally, optimal paths can be obtained for small D by evaluating the path integral 
( pT0|) by the steepest descent method. It follows from Eqs. (Q, (TTOj) that the optimal paths 
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provide the minimum to the functional 

Iff 00 « f l f 

Tl[q(t),\(t),f(t)} = - J J dtdt' f( t )T(t-t')f(t') + jf dtX(t)[q-K(q;t)-f(t)}, (11) 

where \(t) = —ik(t); one can think of \(t) as of a Lagrange multiplier that relates to each 
other the optimal realization of the random force and the path of the system. 

Variational equations for the trajectories that provide an extremum to the functional (fTTT) 
are of the form 

dt'f(t-t')f(t')-\(t) = 0, X(t) + K\q;t)X(t) = 0, q(t) - K(q; t) - f(t) = 0. (12) 

In the problem of the stationary probability density for the system to be in the state qf 
at the time tf (this probability density is periodic in tf with the period T), the boundary 



conditions for Eqs.( fT2|) take the form (cf. 191 ]) 

f{t) for t -> ±oo, A(t) ^ for t -> -oo, A(i) =0 for i > if, 



q(t)^q {0) (t) for i- -oo, q(t f ) = q f . (13) 



In deriving the boundary conditions ( U3j) for t — > — oo we took into account that the system 
fluctuates about the stable state for a long time before the large fluctuation starts (cf. 
Fig. 1). Respectively, one may set in (Q ti — > — oo, g(tj) = q^°\ti) (this is consistent with 
Eqs. (1121) . On the other hand, the motion of the system after it has reached the point qf is 
not important for the large fluctuation. Therefore the constraint on f(t) is lifted for t > tf. 
Clearly, the force decays to zero for t > tf. 

The boundary conditions should be modified if one considers the problem of escape from 
a metastable state. By generalizing the arguments [3] to the case of a non-polynomial 
reciprocal power spectrum l/<3>(u;), one can show that the optimal escape path corresponds 
to the system approaching the unstable periodic state q u (t) (the saddle state) for tf — > oo, 
and X(tf) — > for tf — > oo, because in this case J Q T dt K'(q u (t);t) > 0. The analysis of the 
escape problem is beyond the scope of the present paper. 

Eqs. (1121) . ( 1T31) provide a complete set of equations for the interconnected optimal fluc- 
tuational paths of the system and the force, qopt(t\Qf,tf) an d fopt(t\qf,tf), for reaching the 
state qf at the instant tf. The probability to reach this state, according to Eq. ffTUl) . is of 
the form 



p(q f ,t f ) oc exp[-R(q f ,tf)/D] , R(q/,tf) =K[q opU \ opt J opt } = minfc[g, A, f). (14) 
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where X op t(t) is the optimal Lagrange multiplier as given by Eqs. (fT2l) . (|T3|) . 

We note that Eqs. ffl2l) . (TT3l) may have several solutions. In this case, the physically 
meaningful solution is the one that provides the absolute minimum to the functional 1Z. The 
criterion of applicability of the approach is R/D 3> 1 - it is this condition that determines 
how small should the noise intensity D be. 

An important feature of fluctuations induced by non- white noise, as it is clear from (|T2l) . 
ffT3]) . is that the optimal force f op t(t) does not become equal to zero once the system has 
reached the state qf. Time evolution of the optimal force is given by the equation 



fopt{t) = / dt'w - f)\ opt (t'), 0(t) = zrV(t). (is) 

J — oo 

The function (p(t) is the noise correlation function rescaled so that it was independent of the 
characteristic noise intensity D (the noise correlator oc D, cf. ([I]), ([5])). 

Based on Eq. f|T5|) one can predict how the system will move, most likely, after the state 
qf is reached. The trajectory of the system is described just by the equation of motion 
with the random force given by Eq. ( 1T51) . We note that, even for the state qf lying in the 
basin of attraction to the initially occupied stable state from which the fluctuation starts, 
on the way back from qf the system may come not necessarily to the same stable state, but 
to a different state. This is in contrast with what happens in systems driven by white noise, 
unless qf lies very close to the basin boundary. An example of such fluctuations in systems 
driven by colored noise was, in fact, considered in [ljj, although trajectories of the system 
after it has fluctuated to a remote state were not analyzed. 



C. Vicinity of the stable state 

To illustrate the solution of the variational problem (fT2|) . f[T51) we will consider a simple 
case where the state qf is close to the stable state q^(tf), so that the force K(q;t) can 
be linearized in q — q^°'(t), and yet the difference \qf — q^(tf)\ is big enough so that the 
asymptotic expression (fT^j) applies. With account taken of Eq. ([7j), we obtain: 

\(t) = u{t,t f )\ f , q(t)-q<-°\t)= f drf(r)u(r,t), (16) 

J — oo 

I drK'{q^{T);r) . 



u(t, t') = exp 
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One can see from Eq. (fT6l) that 

\f = g- 1 (t f ,t f )[q f -q« s >{t f )], g(t,t')= f dr f dr' 4>{r - r')u(T,t)u(r f ,t'). (17) 

J — oo J — oo 

From Eqs. (|15p - (|17|) one can find the activation energy f[T^|) of reaching the state qf, 

R(qf,t f )= 1 -g-\t f ,t f )[q f -q^(t f )} 2 . (18) 

The activation energy (fT8|) is quadratic in the distance between the state qf and the stable 
state q(°'(tf). The proportionality factor g~ 1 (tf,tf) depends on the shape of the noise power 
spectrum, and also on the dynamics of the system in the absence of noise. We note that, 
if the regular force K(q; t) is independent of t, the stable state q^ is independent of t, and 
the function g(tf,tf) is independent of tf, as expected. 



III. THE PREHISTORY PROBLEM 



The distribution of paths for large fluctuations can be investigated and visualized through 
the analysis of the prehistory probability density, Ph(qh, th \ Qf,tf) H- This is the conditional 
probability density for a system that (i) had been fluctuating about q^ (t) for a time greatly 
exceeding the relaxation time of the system t re \ and the correlation time of the noise t corr , 
and (ii) arrived to the point qf at the instant tf, to have passed through (and been observed 
at) the point qh at the instant th, th < tf. Using the path- integral formulation ( jTUi) one can 
write the prehistory probability density as 



Ph{Qh,th\ Qf,tf) 



M J Vf(t)V[f(t)} J V^- p 'Vq(t)6[q(t 
Mt) 



Qh] 



(19) 



x exp 



dt 



[q-K{q-t)-f{t)]--K\q-t) 



dqph{q,t\q f ,t f ) = 1. (20) 

The normalization constant M in the general expression for ph f|T9|) is defined by the condi- 
tion (l20l) . Throughout this section we assume that tj — > — oo. 

We expect (and will confirm a posteriori) that, for small D, the distribution ph{q, t\ qf,tf) 
peaks sharply for q lying close to the optimal fluctuational path q op t(t\qf,tf). Therefore in 
evaluating ph one can expand the exponent in V[f] and the term with k(t) in the exponent 
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in Eq. ([19]) in the deviations 5f(t),Sq(t),Sk(t) from the optimal realizations of the random 
force f op t{t), the trajectory q op t(t), and k opt (t) = iX op t(t)- 

In this paper we will not address the problem of singularities of optimal paths in systems 



driven by colored noise [19|, which has been understood only recently for white- noise driven 
systems If the point (qf, tf) is far from singularities, it suffices to keep in the aforemen- 
tioned expansion only the second-order terms in 5f(t),5q(t),5k(t). Integrating over 5f(t) 
and writing S[q(th) — qh] as an integral, one can re-write Eq. ffl9|) in the following form 

Ph (q h ,t h \q f ,t f )=M 1 J°° j V^- J Vq(t) exp (S[6k(t), 6q(t)]/D) , (21) 

where the quadratic functional S is given by the expression 

S[5k,5q] = \ 11 1 dtdt' 6k(t)j>{t - t')5k(t') 

2 J J-OO 

V 1 /"*/ 



-i j dt5k(t)[5q(t) - K'(t)5q(t)] - - j dt\ opt (t)K"(t)5q 2 (t) + ia[q(t h ) - q h }(22) 



with 



5k(t) = k(t) - k opt {t) = k(t) - iXop t (t), 5q(t) = q(t) - q op t{t), 
K'{t) = K'( qopt (t);t), K"{t) ee K"{ qopt {t)-t) = d 2 K/dq\ 

and with the boundary conditions 

5q{-co) = Sq(t f ) = 0, Sq(t h ) = q h - q opt (t h ), 6k(—oo) = 0. (23) 

It is convenient to rewrite the expression for S 1 in the matrix form, 

s , HVi,^] = J /"/"' dtd? (Mt),Mt))n(t,t') | fia[?(t h ) -?/,], 



2 



where Ti is a Hermitian operator, 



cS(f - f) (A"(C) + d/M ) -*(( - t')\^{t)K"(t) 

In obtaining Eq. fT25|) we took into account the boundary conditions (T231) . As we will see, 
of immediate interest is the value of S[5k, Sq] for purely imaginary 5k and real 5q, which 
justifies the unusual form of S 
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A. The variance of the prehistory probability distribution 



Integration over 5k(t),5q(t) in the expression (121]) can be performed by the steep- 
est descent method. In this method, one has to find the extremum of the quadratic 
functional S, which requires solving the following equations for the extreme values of 
Sk(t) = -i^ 1 {t),8q{t)=^ 2 {t): 

dt"H(t,t') I ^ ) ) = -ia8(t h \, Mth) = Qh ~ q op t(t h ) (26) 



M?) J V 1 

where the value of a is determined by the boundary conditions (123]) . 

The solution of Eq. (126]) can be sought in terms of the Green function Gij(t,t') which 
provides the solution to the equation 

/ 1 dt'H vll (t,t')G ll3 (t',t") = 8 i5 8{t - t"), Gij&t') = Gji(lf,t), G a (t,t f ) = 0. (27) 

J — oo 

Here, summation is performed over repeated subscripts i±; the subscripts take on the 

values 1,2. The functions ipi in Eq. (126]) are expressed in terms of the function G as 

^•(*) = [Qh ~ q op t(t h )]G j2 (t, t h )/G 22 {t h , t h ). (28) 

Clearly, the functions ipj are proportional to the distance of the point qh, for which the 
prehistory probability density is sought, to the optimal fluctuational path q op t{th\qf,tf). 

Using Eq. (I24p . the prehistory probability density can be expressed in terms of the func- 
tions ^1,2, and then in terms of the Green function G, 

a/4- \„ + w-in ( \qh - qo P t(th\qf, t f )] 2 



p h (q h ,t h \q f ,t f ) = (2irD<rtft h \q f ,t f ))- l ' 2 exp [- ^ Z) ' , V (29) 

vKthlqf^f) = G 22 {t h ,t h ). 

The distribution ph is Gaussian, with a maximum on the optimal path. It is seen from 
Eqs. ( 128]) . ( 129]) that the variance of the distribution is given immediately by the component 
of the Green function G 22 (th,th). For very large tf — th the variance becomes independent 
of tf, and Eq. (129]) describes the stationary distribution about the stable state. 

Alternatively, the problem of the variance of the prehistory probability distribution can 
be solved in terms of the eigenfunctions and eigenvectors of the appropriate Hamiltonian. 
For white-noise driven systems this method was discussed earlier 
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We note that the above analysis makes it also possible to investigate the "post-history" 
probability distribution: the distribution of the paths of the system after it has fluctuated 
to a remote state. For white-noise driven systems, this distribution peaks at the path along 
which the system, prepared initially in the state qf, goes down to the stable state q(°\ as 
clearly demonstrated experimentally by Luchinsky and McClintock joj. As discussed below 
Eq. (fT5|) . colored noise is not turned off once it has driven the system to a given state, and 
therefore it affects the motion of the system after the state qf has been reached. 

Formulation of the post-history problem requires changing in all above equations to inte- 
gration over the paths q(t) and the auxiliary field k(t) for t varying from — oo to oo (instead 
of — oo to tf), with the condition that the paths go through the state qf at the instant tf. 
The final answer is again given by Eq. f[2T?j) . and the distribution peaks at the most probable 
path for the motion of the system after the state qf has been reached. 

IV. PREHISTORY PROBABILITY DISTRIBUTION CLOSE TO THE STABLE 
STATE 

Explicit expressions for the prehistory probability density can be obtained if the final 
point qf lies close enough to the stable state q(°'(tf), in which case the equations for optimal 
paths of the system and the force are linear. In the prehistory problem, to the lowest order in 
\qf — q^(tf)\, one can neglect the term in (|25|) with X opt K" oc Qf — Q^{tf) (cf. ffT7|) ). and can 
also replace K' by its value for the stable state. Eq. (1261) may then be immediately integrated. 
After straightforward but somewhat tedious calculations one obtains the following expression 
for the reduced variance of the distribution p h : 

4(th\q f , t f ) = [g(t f , t f )g(t h , t h ) - g 2 (t f , t h )] /g(t f , t f ) (30) 

(the function g(t,t') is defined in Eq. (ITT])). 

We emphasize that Eq. (I3"0j) applies for an arbitrary shape of the power spectrum of the 
noise (however, we assume that the function 0(r) decays at least exponentially for |r| — > oo). 
It applies also for an arbitrary periodic driving. It may be further simplified in the absence 
of periodic driving, in which case K' = —a = const, with a > 0, and we obtain from (ITT)) 

g (t,t') = — [ diuD^(uj)(a 2 + u; 2 )~ 1 e Mt ~ t ' ) - (31) 
2tt J 
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This expression makes it simple to calculate the variance of the prehistory probability density 
for an arbitrary shape of the noise power spectrum $(u;). The results of these calculations 
for noise with the Gaussian power spectrum centered at zero frequency are shown in Figs. 2. 
It follows from this figure that the noise color changes the broadening of the prehistory 
distribution very substantially. An important consequence of Eqs. (1301) . (1311) is that the 
variance a 2 (th\qf,tf) oc (£/ — th) 2 for small t/ — th- This is qualitatively different from the 
linear time dependence of a\ in the case of white noise, known from 7, sj]. The reason is 
that, for systems driven by colored noise, the mean square displacement over the time At, 
which is small compared to the noise correlation time t corr , is proportional to (At) 2 , not to 
At, as for white-noise induced diffusion. 




FIG. 2: Time dependence of the scaled variance a 2 of the prehistory probability density in a linear 
system driven by noise with the Gaussian power spectrum $(w) = D exp(— u; 2 /2a 2 <7 2 ), for different 
values of the dimensionless variance of the noise spectrum a 2 . The time is scaled by the decrement 
of the system a = —K'(q^). 



For large tj—th, which exceeds the relaxation time of the system and the correlation time 
of the noise t re i,t corr , the prehistory probability distribution goes over into the stationary 
distribution which is described by the function R(qh, th) = [qh~q^ (th)] 2 /2er^, where R(qh, t h ) 
is given by Eq. ffTS]) . 
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V. CONCLUSIONS 



In contrast to white noise driven systems, for colored noise, after the noise has driven 
the system to a remote state g/, it does not become small at a time. As the noise decays it 
drives the system further along a certain path. This path differs from the path which the 
system would follow if it were prepared in the state qf "by hand", not as a result of the 
fluctuation. 

In the present paper we provided a general formulation of the problem of large fluctuations 
in systems driven by Gaussian noise. This formulation makes it possible to describe optimal 
fiuctuational paths of the system, and also to evaluate the width of the tube of fluctuational 
paths that arrive at a given target state. The latter is done using the prehistory probability 
distribution for the system to have passed through a given point on its way to the target 
state. The tube of the fluctuational paths is centered at the optimal path. Evaluation of the 
width of the tube has been reduced to solution of a linear equation. Explicit results have 
been obtained for the fluctuations in the linear range close to the stable state, and the effect 
of noise color in this domain has been analyzed. 
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